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Abstract. We calculate the three-point correlation function evaluated at horizon 
crossing for a set of interacting scalar fields coupled to gravity during inflation. This 
provides the initial condition for the three-point function of the curvature perturbation 
in the Sasaki-Stewart 6N formulation. We find that the effect is small, of the order of 
a slow-roll parameter, and that the non-gaussianity can be determined on large scales 
once the unperturbed background evolution is known. As an example of the use of our 
formalism, we calculate the primordial non-gaussianity arising in a model of assisted 
inflation. 
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1. Introduction 

The inflationary paradigm is an attractive proposal for the evolution of the very early 
universe [U El El El E] ■ During inflation the universe underwent a phase of accelerated 
expansion driven by one or more self-interacting scalar fields. Despite the appealing 
simplicity behind the central idea of inflation, it has proved difficult to discriminate 
between the large number of different models that have been developed to date jSj. 
The simplest classes of models typically predict that the spectrum of primordial density 
perturbations generated quantum-mechanically during inflation should be nearly scale- 
invariant and Gaussian-distributed. (For a review, see, e.g., Ref. [Jj.) These generic 
predictions are in good agreement with cosmological observations, most notably those 
arising from measurements of the temperature anisotropy and polarization of the Cosmic 
Microwave Background (CMB) [HIE]. Such observations lead to strong constraints on 
theoretical model building. 

However, if further progress is to be made in our understanding of the microphysics 
of the early universe, it will become increasingly necessary to extend the theoretical 
framework beyond the leading-order effects of scale- invariant, Gaussian fluctuations. In 
particular, deviations away from Gaussian statistics represent a potentially powerful 
discriminant between competing inflationary models and have attracted considerable 
recent interest [TQl [TTl H21 USl HSl HEl HZl HH1 [THl - Moreover, measurements of CMB 
anisotropies are now approaching a level of precision where such differences could soon 
be detectable and, indeed, the first non-gaussian signals in the data may already have 
been observed [HI EH E2 EH| (see also gl). 

In general, any Gaussian random variable X with zero mean and variance a 2 has a 
probability measure given by 



The rules of Gaussian integration imply that (X n ) vanishes whenever n is odd, whereas 
(X n ) can be written in terms of (X 2 ) for even n. Accordingly, we expect violations of 
Gaussian statistics to manifest themselves as deviations from these simple rules. The 
same principle holds when dealing with fluctuations Sip in some quantum field, where 
the semiclassical probability density is given by ¥(5(f) ~ e _/ ^ and I is the Euclidean 
action. When / describes a free theory with no interactions, 5(p may be viewed as a 
Gaussian random variable. However, if interactions are present, deviations from pure 
Gaussian statistics will arise. More specifically, if 8ip represents an 'almost Gaussian' 
random variable, then it is expected (in a largely model-independent sense) that the 
three-point function (5cp 3 ) will provide the largest contribution to the non-gaussianity. 

In typical single- or multi-field inflationary models, the level of non-gaussianity 
can be calculated analytically [10]. The single- field scenario has been well-studied 




(1) 



and the correlation functions of X have the form 




(2) 
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[251 1213 HH Uni ITU] , and some effort has focused recently on non-standard theories 
of gravity, in the hope of identifying observational signatures that differ from those 
of canonical inflation based on the Einstein action [211 EBl CHI E|- Some earlier 
investigations considered non-gaussian fluctuations, but were specifically directed at the 
formation of primordial black holes during the final stages of inflation, when the slow- 
roll approximation has broken down [23J EDI- I n general, however, the superhorizon 
evolution of the curvature perturbation during single-field inflation is simple, and there 
is a robust prediction that the level of primordial non-gaussianity as measured by the 
three-point correlation function cannot be large [Tl] , even allowing for exotic effects such 
as higher-derivative operators of the scalar field or a speed of sound that is different from 
that of light [HE]- 

After the end of inflation, one must evolve this small non-gaussianity through 
the epochs of radiation and matter domination in order to arrive at the temperature 
anisotropy 5T/T which is seen on the microwave background sky pH] 13*2*] . This 
evolution after the end of inflation gives an additional source of non-gaussianity which 
is comparable to the primordial non-gaussianity in the case of single-field inflation. 

By comparison, our understanding of the multi-field scenario is less developed, 
although some progress has been made [2niE2ll3ll3liaEEl[IBl[^ 
Typically, the evolution of the curvature perturbation on superhorizon scales implies that 
a significant non-gaussian signal can in principle be generated during inflation. Indeed, 
a fairly general prescription for calculating the three-point function of fluctuations was 
recently presented ^H] , in which it was assumed that the curvature perturbation evolves 
substantially outside the horizon due to pressure effects after the end of inflation. 
Such superhorizon evolution effects generically arise when more than one fluid is 
cosmo logically important, or when a single fluid is present but with an indefinite 
equation of state. In the multiple-field case, therefore, the primordial three-point 
function consists of a microphysical component, calculated from the properties of the 
matter fluctuations around the epoch of horizon exit, and a superhorizon component, 
which is determined purely by large-scale gravitational effects. If the superhorizon 
component becomes large, which typically will not happen until after the end of inflation, 
it may dominate the primordial non-gaussianity, in which case it is consistent to neglect 
the microphysical contribution. This actually occurs, for example, in the curvaton 
scenario [J21 El EJ El El- I n this respect, the theory has some similarities with the 
stochastic approach to inflation, which has been revived recently within the context of 
non-gaussianity [THl E] • 

However, in the absence of any estimate for the microphysical contribution to the 
three-point function, the above prescription is necessarily incomplete, since it is not 
possible to verify explicitly that the superhorizon contribution dominates. In this paper, 
we supply the missing ingredient by calculating the three-point function evaluated at 
the epoch of horizon crossing for a generic set of M scalar fields {tp 1 } coupled to gravity. 
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We consider a (Lorentzian) action of the form 

S= l -Jd A x^g~[R- Qu d^d^ J - 2V{<p)] . (3) 

where the interaction potential V(ip) is assumed to be an arbitrary function of the 
fields, ip 1 , and Qu represents the metric on the manifold parametrized by the scalar 
field values. (We refer to this metric as the 'target space' metric). The action for 
M canonically normalized scalar fields minimally coupled to gravity is recovered by 
specifying Qu = 8u, where Su represents the Kronecker-delta. More generally, the 
kinetic sector of (JHJ) may be invariant under a global symmetry group J, in which 
case the scalar fields parametrize the coset space J/K with a line element given by 
ds 2 = Gu dp> I d(p J , where K is the maximal compact subgroup of J. In such models, 
the metric components generically exhibit a direct dependence on the scalar fields (p 1 . In 
this paper, we restrict our attention to the case where the target space is approximately 
flat, at least over the region dynamically explored by the fields ip 1 , and the metric 
can be brought to the field-independent form Qu = Sij by an appropriate choice of 
parametrization. On the other hand, we allow the potential V(<p) to be quite arbitrary, 
subject to the usual slow-roll conditions. This choice of action covers a wide class of 
multi-field inflationary scenarios. 

The outline of the paper is as follows. In Section El we briefly describe how 
to calculate the primordial non-gaussianity produced during multi-field inflation, as 
measured by the three-point function of the curvature perturbation. The crucial 
ingredients for the calculation are: (i) an expression for the curvature perturbation 
itself, which is provided by the Sasaki-Stewart 5N formalism and: (ii) an estimate 
of the scalar three-point function just after horizon exit. This three-point function is 
derived in Section El and is presented in Eqs. (|58 j) -(p§ J) . It represents the central result 
of the paper. In Section we verify that our formalism reproduces the well-known 
result of Maldacena jTI] when specialised to the case of a single scalar field. In Section |3] 
the non-gaussianity generated during assisted inflation E3 HE] is calculated as an 
example of our method. Finally, we conclude in Section |3 

Units are chosen such that c = h = Mp = 1, where Mp = {&nG)~ l l 2 represents the 
reduced Planck mass. 



2. Non-gaussianity from multiple fields 

2.1. The background model 

We assume throughout that the background model corresponds to the spatially flat 
Friedmann-Robertson- Walker (FRW) universe, with a metric 

ds 2 = -dt 2 + a 2 5 ij dx l dx ] . (4) 

The scalar field equations take the form of a set of coupled Klein-Gordon equations 

HHJUEEHJEI!: 

(p 1 + ZHtp 1 + Q IJ Vj = 0, (5) 
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where H = a/ a is the Hubble parameter, a is the scale factor, ip 1 = dip 1 /dt and a 
comma denotes partial differentiation. 

The gravitational equation of motion is 

2H + 3H 2 = ~Gutf<p J + V(<p), (6) 
together with the Friedmann constraint equation 

3H 2 = ^g IJ <p I <p J + V(ip). (7) 
Eqs. (0) and (JJJ) together imply that 

H = -iftj^^. (8) 
There is also an analogue of the single-field Hamilton- Jacobi equation, which reads 

Hi = ~<pi. (9) 

As in single-field inflation, it will prove convenient to introduce a set of slow-roll 
parameters. One such parameter is the slow-roll matrix e IJ : 

TJ < f I< f J nr IM r JN H M H ,N T J nn x 

" = JlP =2g 9 ~ll^ =££ > (10) 

where we have used ©, and e 1 is the vector 
V2H H 

The standard slow-roll parameter, e = —H/H 2 , can then be expressed in terms of e IJ , 
since e = tre IJ = QijS ij ' . As a result, we generally expect \s IJ \ <C 1 whenever slow-roll 
is valid except in models where finely-tuned cancellations between the components of 
e IJ may occur. Specifically, for generic theories which are not tuned to have anomalous 
parameter values, one should expect 

*" = <>(£), *' = o(^). (12) 

The time derivative of e IJ is determined by 

e IJ = 2eH (e IJ - r] IJ ) , (13) 



where we have defined a second slow-roll matrix 

v u = pV + pV (14) 

AHH 

This matrix generalizes the slow-roll parameter 7] of single-field inflation, as can be 
deduced by substituting in Eq. (jHJ) to yield if* = —ip/Hip = rj. 
Finally, we introduce a third matrix defined by 

Vu = (15) 

This is related to e IJ and 7] IJ by the simple rule 

nM{I f J)n 

£ " + v " = fj MN y (16) 



Primordial non-gaussianities from multiple-field inflation 



6 



The parameter fju represents a generalization of the slow-roll parameter V"/V of single- 
field inflation, and can be related to the alternative definition —<p/H(p. In multi-field 
models, however, the index arrangement in (jlfij) means that it is convenient to maintain 
two separate definitions. 

We note that Eq. (jlfij) implies that rj IJ and fjjj are of the same order in slow-roll. 
The utility of this formulation of the slow-roll approximation is essentially restricted to 
the case of a flat target space. 

It will also be useful to have on hand an explicit expression for the number of e- 
foldings in this model. One defines the number of e-folds N which elapse between some 
initial value cij of the scale factor and a final value a/ as iV = lna//aj, so 

diV = Hdt = — dlntf = — %d(p T . (17) 
2.2. The uniform density curvature perturbation 

We now consider perturbation theory around the above homogeneous background. Since 
it is assumed that the scalar fields (p 1 dominate the energy density of the universe, 
any perturbations dip 1 in these fields necessarily produce a disturbance in the energy- 
momentum tensor, which back-reacts on the spacetime curvature. Therefore, when 
working with the perturbation theory of these scalar fields, we must also take into 
account the effect of scalar metric fluctuations in order to include the effect of this 
back-reaction. 

To first-order in perturbation theory, these metric fluctuations can be written as 



ds 2 = -(1 + 2$)dt 2 + 2a 2 Ei dxMt + a 2 [(1 - 2^)<% + 2-E I y] dx l dx j . (18) 

Not all of these degrees of freedom are dynamically important, however, since they are 
related by the gauge transformations t i— > t + 5t and x % i— > x l + Sx l , which reshuffle $, 
B i and E^j. 

To obtain a gauge-invariant measure of the strength of such fluctuations, it is 
sufficient to consider scalar curvature invariants of the three-dimensional hypersurfaces 
t = constant. The simplest such invariant is given by 

y=^ d -2 R (3) y (19) 

where is the three-dimensional Ricci scalar on the spatial hypersurfaces, and 
d~ 2 is the inverse Laplacian^. In single-field models, this quantity has an attractive 
interpretation. Under a temporal gauge transformation t i— > t + H5t, the field \l/ 
transforms according to the rule 

^ ^ ^ + H5t. (20) 

| This is most simply defined in Fourier space, where the Laplacian d 2 is represented as multiplication 
by —k 2 . Similarly, the inverse Laplacian corresponds to multiplication by —k~ 2 . 
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Hence, the curvature scalar on comoving hypersurfaces (defined by the condition that 
the scalar field perturbation 5<p vanishes) is given by 

n = v 5v=0 = y + H 6 -^. (21) 

One may also consider the curvature scalar on uniform density hypersurfaces, 
defined by 5p = 0. This scalar is denoted by (, where \(\ = |^5p=o| jSH], but there 
is no consistent agreement on sign conventions. For adiabatic perturbations generated 
during single- field inflation, the two sets of hypersurfaces coincide [HB] , so modulo choices 
for signs, it follows that TZ — \(\. One can freely work in either the comoving or uniform 
density gauges, and the gauge-invariant variable TZ — |£| = ty+HSip/ip is the only degree 
of freedom in the theory. This variable mixes the fluctuations arising from the metric 
and the scalar field, neither of which has a separate gauge-invariant interpretation. 
Consequently, at the quantum level one considers the creation of TZ- or (^-particles from 
the vacuum rather than that of inflaton particles or gravitons. 

After horizon exit, TZ asymptotes towards a time-independent quantity and may be 
viewed as a function only of position. It can therefore be absorbed by an appropriate 
choice of coordinates on the spatial hypersurfaces. This implies that distant regions 
of the universe differ from each other only by how much one has expanded relative 
to the other. Since this is a symmetry of the equations of motion, TZ is conserved on 
superhorizon scales [TTj. 

In multi-field models, the situation is not so simple. It is still convenient to 
work with the uniform density curvature perturbation, (, which on large scales is 
still equivalent to the comoving curvature perturbation TZ. However, expressing ( as 
a mixture of the metric perturbation \l/ and the scalar field perturbations Sep 1 is much 
more involved. This was first achieved at a linear level by Sasaki & Stewart [15]; who 
used fl2Hj) to express the difference between an initial uniform curvature hypersurface 
and a final comoving hypersurface as 

(=\Hdt\=dN, (22) 

where N is the integrated number of e-folds, as defined in ([TYf) . On sufficiently large 
scales we expect (given suitable assumptions about the dynamical behaviour which 
permit us to ignore time derivatives of the perturbations) that each horizon volume will 
evolve as if it were a self-contained universe, and it therefore follows that (as shown 
beyond linear order by Lyth, Malik & Sasaki [53]) 

C = Nj V + ^N u V 6<p J + ■■■, (23) 

where the dip 1 express the deviations of the fields from their unperturbed values in some 
given region of the universe, and we employ the standard summation convention for 
the scalar field indices J, . . . , J. It is this expression for ( which was used by Lyth & 
Rodriguez to express the primordial non-gaussianity of ( in multi- field models [Ej. (See 
also [53].) 
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When more than one field is dynamically important, Eq. (J23|) allows TZ to evolve 
outside the horizon while inflation is still occurring§, in contrast with the single-field 
case. It can be shown that the evolution of TZ is sourced by a non-adiabatic pressure 
perturbation 

TZ ~ — <fy nad , (24) 

P + P 

where £p n ad can be written in the form 

s Pmi=i> r^p( 6 -E- S 4), (25) 
V v p J 

and T represents the so-called isocurvature perturbation that measures the displacement 
between hypersurfaces of constant density and of constant pressure. (Note that 
T vanishes whenever p is a definite function of p.) Thus, whenever isocurvature 
perturbations are present, which are essentially inevitable in a general multi-field 
scenario, TZ will evolve outside the horizon while inflation is still occuring, and this 
evolution must be taken into account in the inflationary prediction I50j . The Sasaki- 
Stewart formalism conveniently handles these complexities for us. 



2.3. The three-point function in multiple-field models 

In the remainder of this section, we briefly review the results of [TS] . 

In order to connect with observations, any prediction for the non-gaussianity 
measured by the three-point function must be expressed in terms of an experimentally 
relevant parameter. A common choice is the non-linearity parameter, /nl ; which 
expresses the departure of £ from a Gaussian random variable: 

C = Cg ~ jj/NL * & (26) 

where ( g is Gaussian and we have written a convolution product since /nl is typically 
momentum dependent (and defined by multiplication in Fourier space [SZ1 EH])- We 
have dropped an uninteresting zero mode in ([26)1 . which only serves to fix (£) = 0. 
This so-called primordial non-linearity is processed into an observable non-linearity of 
temperature anisotropies /^l as described in 

To relate /nl to the three-point function, we follow Maldacena JT] and write the 
latter as 

47T 4 

(C(k 1 )C(k 2 )C(k 3 )) = (2nf5(k l + k 2 + k,)-—\Al\ 2 A, (27) 

Hi K i 

§ To be clear, we should point out that there is evolution of the metric perturbations after the end 
of inflation. One source of evolution, arising from the evolution of perturbations during the epochs of 
radiation and matter domination which follow reheating \\\\ , is universal to all models of inflation, 
irrespective of whether many fields are important or only one. Additionally, in multiple field models, 
N will typically continue to vary after the end of inflation. Since (|2,'iH is still valid, this means that £ 
will continue to change. In this paper, we only evaluate £ up to the end of inflation. 
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where A is a momentum-dependent function that can be calculated explicitly, and A^ 
is the dimensionless power spectrum of (, defined in the standard way by 

k 3 

(C(ki)C(k 2 )> = (2vr) 3 5(k 1 + ka)P C (*i). A K k ) = ^~2 P dk), (28) 

such that ((() = J °° A^(fc) din k. The non-linearity parameter can now be written in 
the succinct form [TT] ||: 

! - = (29) 

It only remains to use to relate /nl and the scalar field three-point functions. 
We show in Section |3] that the two-point function for the scalars satisfies 

(V(ki)<V(k 2 )> = (27r) 3 5(k 1 + k 2 )^AlQ IJ , (30) 
and that the three-point function can be written in the form 

47T 4 

(V(ki) 5<p J (k 2 )6ip K (k 3 )) = (2n) 3 5(k 1 + k 2 + ks)^ \A 2 fA IJK , (31) 

where A^ is the spectrum of a massless scalar field in de Sitter space and A IJK is a 
momentum-dependent function given in Eq. (J69j) . The determination of A IJK represents 
the principal result of this paper. 

Eq. (J23)) implies that the power spectrum for ( is given by 03] 

Aj5 = AlNjNjG IJ . (32) 

On the other hand, the three-point function which follows from (|23J) can be written as 

(C(k 1 )C(k 2 )C(k 3 )) = iV j/ iV iJ iV^(V(ki)^ J (k 2 )^(k3)) 

+ ^N I N K N MN (8if I (k 1 )8ip K (k 2 )[8^ M * 5<p N ](k 3 )) + perms + • • • , (33) 

where * denotes the convolution product, and for brevity we have omitted the remaining 
cross terms in the expansion of ( 3 . (In particular, there will be terms involving the 
product of a spectrum and a bispectrum, of the form 

^iVziV^iVM^ V(ki)[V" * 5<p L ](k 2 )[5<p M * Sp N ](k 3 )) + perms, (34) 
and a term involving the product of the spectra, of the form 

InjjN kl N mn ([5^ * V](ki)[V * 6<p L ]{k 2 )[6<p M * 5ip N ]{k 3 )). (35) 



The first of these was ignored by Lyth & Rodriguez [T3] , because it is proportional to 
the three-point function for the Sip 1 , which these authors assumed to vanish. These 
terms can all be treated similarly jSSl UH EH| when included in / NL , but some of the 

|| There are other relations between /nl and the three-point function which are common in the 
literature, usually written in terms of the bispectrum, but they are all equivalent to this one. In 
the present case, the use of the dimensionless power spectrum is especially convenient because it will 
allow us to factorize the index structure in the multiple-field three-point functions. 
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intermediate expressions can become rather long. For these reason we suppress these 
higher contributions.) 

The first set of terms in Eq. (J23)) is just the connected scalar three-point function, 
since we have assumed that all tadpole contributions vanish. This is mandatory, and 
equivalent to perturbative stability of the vacuum. Under the same assumptions, the 
second group of terms can be reduced to products of two-point functions. Consequently, 
by employing (|30 ]) -(|3*T ]) . Eq. can be expressed in the form 



47T 4 

(C(k 1 )C(k 2 )C(k 3 )) = (2 7 r) 3 5(k 1 + k 2 + k 3 ) r — 3 \Al\ 2 A (36) 
where Aq is defined by 



A c = A ijk N ;I N i jN ! k + g IM g KN NjN K N MN kl (37) 

i 

Our expression for /nl; Eq. (j2Hl), now implies that the non-linearity in this model, to a 
good approximation, is given by 

5 A ijk NjNjN k 5 g IM g KN NjN K N MN 

~ 6(iVVje jJ ) 2 EA 3 6 {NiNjguy + '"' 1 j 

where by '• • •' we denote the remaining cross terms from which we did not display. 
In principle there will also be other terms coming from the Taylor expansion of ( (in 
Eq. (J23J) ) beyond second order, but we assume that the series converges sufficiently fast 
that the effect from these higher-order terms is small. Nevertheless it should be borne in 
mind that they are present in principle, and in certain models their contribution could 
be significant. In order to express /nl as in (JSHJ), we have exchanged for A;? in (|3l)]): 
this has the effect of introducing the terms involving NjNjG in the denominators of 
Eq. (|38jl . We note that the elements of this expression were first given in [TH] . 

Generally speaking, we expect both displayed terms in (|3*8Jl to be of roughly 
comparable importance at the end of inflation, although examples of models where one 
term is dominant do exist. For example, assisted inflation represents a model where the 
second term vanishes, so /nl arises entirely from the first term. On the other hand, the 
authors of Ref. ^H] focused on the case where the second term is dominant. Note that 
the first term and any higher cross-terms containing the bispectrum A IJK are strongly 
momentum-dependent, whereas the second term and any higher cross-terms which do 
not include the bispectrum are momentum- independent, up to mild scale-dependences 
which arise from renormalization effects lf)9j. 

We proceed in the following section to calculate the three-point function for the 
multi-field action (JSj). 

3. The scalar three-point function 

3.1. Perturbations in the uniform curvature gauge 

The first step in determining perturbations in multi-field inflation is to select a gauge. 
The discussion of the previous section, where was used to express the curvature 
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perturbation (, implies that the relevant gauge should be the uniform curvature slicing, 
where coordinate time t is chosen in such a way that slices of constant t have zero Ricci 
curvature. Such a metric can be written in the ADM form 

ds 2 = -N 2 dt 2 + hijidx 1 + N l dt)(dx j + N j dt). (39) 

This version of the metric incorporates the scalar perturbations $ and written 
down in (|18j) within the lapse, N, and shift, N l . These latter two quantities appear 
as Lagrange multipliers in the action and for this reason it is considerably easier to 
perform calculations with (J3llj) than with the equivalent expression (|18|). The term 
of ()18|) has been removed by a coordinate redefinition. 

One specifies = a 2 (t)5ij in order to obtain spatially flat slices. The scalar fields 
on these flat hypersurfaces are then written as ip 1 = (Pq + Q 1 , where ip^ are the spatially- 
homogeneous background values of the fields and Q 1 represent small perturbations. In 
these coordinates the action (J3J) becomes 

S = ~\f ^Vh (GubPd i ip l d j ip J - 2V(cp))+^ J N~ l Vh (En&i - E 2 + S^V) , (40) 

where v 1 = ip 1 — NWjtp . 

There are two constraint equations arising from this action. The first follows by 
demanding that the action be stationary under variations of N: 

- QuWd^dw 3 - 2V - ± (EijE* -E 2 + Quv'v 3 ) = 0. (41) 
The second constraint arises from the variation in N l and yields 

^-W) ^^Guv'd^, (42) 

where | denotes the covariant derivative compatible with the spatial metric hij. These 
constraints can be solved to first order by specifying 

N=l + ^GiMQ J , 

Ni = 0"V = ^d- 2 g u (q 1 ^ - ^Q J - §^Q J J • (43) 

It turns out that the higher-order pieces in N and N l are not required, since they cancel 
out of the second- and third-order terms in the action (|4()jl TTllllij . In what follows, we 
will have no need to refer to the total field ip 1 = ipQ + Q 1 , and for notational convenience 
we therefore drop the subscript '0' on (p^ and simply identify ip 1 as the homogeneous 
background field. 

3.2. The second-order theory 

After integrating by parts in the action (J4~U|) and employing the background field 
equations to simplify the result, we obtain the second-order action HH]: 



S 2 = i J dtd 3 x a 3 (^GljQ i Q j - ■^QijdQ I dQ J - MuQ'Q^ 



(44) 



Primordial non-gaussianities from multiple-field inflation 



12 



where a shorthand notation 8Q I dQ J is adopted for the scalar product 5 lJ diQ 1 djQ J , and 
the mass-matrix, Aiu, is defined by 

m " = Vjj - m ■ < 45 » 

Eq. (|44|) is exact and is valid for arbitrary scalar field dynamics. After changing 
to conformal time, r) — J dt/a(t), and introducing the canonical Mukhanov variable 
u 1 = aQ 1 , one obtains a matrix field equation: 



u If > + 



a" 



k 2 - — 

a 



5 J j + M T ja 2 j u J = 0, (46) 



where a prime denotes derivatives with respect to rj. Since the mass term Ai 1 j is not in 
general diagonal, Eq. (j4*U|) does not necessarily factorize into Af independent equations 
for the u 1 . Instead, the interactions described by M. 1 j typically couple the scalar fields 
to one another. 

In the single-field scenario, the mass-matrix Ai 1 j can be expressed as a combination 
of slow-roll parameters. It can therefore be neglected to leading order in slow-roll. To 
arrive at a similar result in the multi-field case, one can express Aiu directly in terms 
of the slow-roll matrices (|10 |l -()16j l . It follows that 

-L.M/J = - 6e u - ee u + erju (47) 

and Eq. (|47|) generalizes the familiar result from single-field inflation. Since we are only 
interested in leading order effects, we can consistently neglect the slow-roll contributions 
arising from Aiu and consider only terms of 0(1). This is the standard approximation 
which is invoked when estimating the amplitude and spectral index of inflationary 
fluctuations to lowest order in slow-roll |HT| I7j. and when calculating the leading-order 
effects of non-gaussianity [TT]. It results in the de Sitter-Mukhanov equation: 

u 1 " + (V ~ ^ V = u 1 " + [k 2 - ^\ u 1 = 0. (48) 

Applying the above simplification implies that we can work with the propagator 

(Q I (x 1 )Q J (x 2 )) = g IJ G*( Xl , x 2 ), (49) 

where G± represents the conventional scalar de Sitter Green's function for a massless 
scalar field [5^] : 

Gt[xuX2) = El x / (1 " **K1 + i^)e-«-») „ > r te 
K ' ; 2k 3 1 (1 — ifc77!)(l + ifc77 2 )e lfc ^ 1 ~ ??2 > m > Vi 

By taking the coincidence limit of the propagator, we obtain the familiar power spectrum 
on large scales (where k — > 0): 

(Q^kOg^ks,)) = (2 7 r) 3 5(k 1 + k 2 )^ 7J ^, (51) 
and it therefore follows that 

A! - (52) 
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The leading-order slow-roll term in the third-order action is given by^f 
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dt d 3 a; a 3 ( - ^Gu&tydQ J - -^Gmn^ M Q N GijQ I Q J 
\ a z AH 

1 1 



a 4 AH 

where ip was defined in 



QiviN i P M Q N Qi.jdQ I dQ J 

and is given to leading order by 



IAJ 



2H 



(53) 



(54) 



In exact analogy with the single-field calculation, it proves most convenient to integrate 
by parts in the term containing dip. One then finds that 



dtd 3 x aQuCfd^dQ 3 = 

dtd 3 x(^ - ^aHG IJ ijdQ I dQ J - ^aGuipdQ'dQ- 7 + aQuQ^&Q^tt) 

Eq. (p>f|) may then be differentiated to determine ip. Keeping only leading-order 
terms, we find that 

ip = -a 2 g u ^d- 2 Q J - ^Gu^d~ 2 Q J . (56) 

An 

It should be emphasized that should not be employed directly in the action 
since it involves second derivatives of the canonical field Q 1 and therefore changes the 
order of the field equations. Instead, by neglecting subleading terms in slow-roll, one 
may use the first-order equation of motion 



1 5L 



a 3 5Q 1 



-3HQ I -Q I + —d 2 Q I + 0(e) 



to eliminate the second- derivative term Q 1 in (|5fij) . It follows that 



2H 



<p 2 d- 



1 5L 



a 3 5Q 1 



3HQj + —d 2 Q! + 0(e 



(57) 



(5f 



Substituting (J58|) into (J53J) and performing a further integration by parts where necessary 
then results in the equivalent third-order action 



dt d 3 x 



^QjQiQ 1 - ^d-'QjQid'Q 1 



+ 



5L 



SQi 



(59) 



The last term in (J59|) . which is proportional to the first-order equation of motion 
SL/SQ^i, is familiar from the single- field calculation, and represents a field-redefinition. 

% This expression should be compared, for example, with Eq. (3.6) of to which it reduces in the 
single-field case. 
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It can be eliminated from the action by transforming the fields Q 1 to new fields Q 1 , 
which satisfy 



where it makes no difference to leading order whether we write Q or Q in the quadratic 
terms. Since we are keeping only terms up to and including third order in Q, this 
redefinition has no effect on the interactions described in (|59jh and we may freely set 
Q 1 = Q 1 there. However, the redefinition does modify the quadratic part of the action, 
Eq. (|44j). which transforms as 



The extra term here exactly suffices to cancel the interactions in (159)) which were 
proportional to the first-order equations of motion, SL/SQ 1 ^. 

Although this field redefinition is extremely convenient for the purposes of 
calculation, we are ultimately interested in the correlators of Q 1 and not Q 1 . The 
correlators are related via the standard prescription 



which is just Wick's theorem. This allows us to add in the contribution of the redefined 
terms in ()60|) at the end of the calculation. 

In the present calculation, the very complicated field redefinition found by 
Maldacena (displayed in Eq. (3.8) of Ref. [TT]) does not appear. This is because such a 
redefinition of fields is almost exactly equivalent to a translation from the comoving to 
the uniform-curvature gauge. In this paper, we have worked in the uniform-curvature 
gauge from the outset, which results in a considerable simplification even for the example 
of a single field. In some sense, this effect demonstrates that the uniform-curvature 
gauge is especially convenient for calculations of this sort: in the single- field case, one 
is free to work equally in the comoving gauge (as in Refs. Ej) or the uniform- 
curvature gauge, but in order to actually perform the functional integrals implicit in 
the three-point correlator (CCC)> one m ust eventually return to the uniform-curvature 
gauge, except for an anomalous field redefinition related to the effect of the Q 1 term. 

It is interesting to note that if we had naively kept Q 1 in the action, instead 
of replacing it with the first-order equation of motion, we would have omitted the 
contribution of the quadratic terms in ([60)1 . This may seem paradoxical, but, as 
we pointed out above, the appearance of Q 1 in the action changes the order of the 
field equations, and consequently their dynamical behaviour. We wish to retain only 
perturbations around the branch of solutions described by (J3j), and no others. For this 
reason, we are obliged to eliminate any second- or higher-derivative terms which appear 
in the slow-roll expansion of the action for the Q 1 , and we achieve this by employing the 



Q 1 = Q! -F I {Q) = Q I 




(60) 




(61) 




(62) 
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background attractor behaviour of the second-order action to select the correct solution 
around which we perturb. 

3.4- Calculating the three-point function 

We are now in a position to calculate the three-point function (Q 1 (ki) Q J (k 2 ) Q K (k 3 )) 
corresponding to the action (|59j) . where we work with the redefined field Q 1 , defined in 
Eq. (ffiOj). so that the 5L/5Qi\i terms are not present. A detailed description of how the 
three-point calculation is carried out has been presented elsewhere [TT1 [TJ] , so here we 
merely record the results. 

• The first interaction to consider is given by 

d V d 3 x f-^QjQ'^A, (63) 

where we have converted the integral to conformal time. This gives a contribution 
to (Q 1 Q J Q K ), after translation to Fourier space, of the form 

- i(2vr) 3 5(^ kl )^^^^g JK J d V k 2 2 kl(l - iAw)e*"> + perms, (64) 

where k t = ki + k 2 + k 3 is the total (scalar) momentum. The integral can be 
evaluated after Wick rotation onto the positive imaginary axis and it then follows 
that this contribution is equivalent to 



ff 3 i 

V 4 EW 



_|3_!tt ) (65) 
_perms ^ * * . 

where an asterisk '*' denotes evaluation at horizon crossing. The sum in this and 
subsequent expressions is over all ways of simultaneously rearranging the indices J, 
J and K and the momenta ki, k 2 and k$, such that the relative positioning of the 
fc's is respected. (In other words, when exchanging indices / and J, for example, 
one should also exchange k\ and k%, and so on.) 

The second relevant interaction is 

/ dT]d3x (-M^ d ' 2Q ' jQ ' ld2Ql ) ' (66) 

where again we have rewritten the integral in terms of conformal time. After 
translating to Fourier space and Wick rotating the resulting integral, we have that 



pp pp 
k t k 2 t 



(67) 

V K t K~t J 

_perms x . 

Once the field redefinition terms in ()59|) have been introduced back into the three- 
point function, using one finds that 

47T 4 

(Q / (k 1 )g J (k 2 )Q^(k 3 )) = (27r) 3 5QTk l ) f — ^\Al\ 2 A IJK (h,k 2 ,h), (68) 

1 Li K i 
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where the spectrum is to be evaluated at horizon crossing, and A IJK is a momentum- 
dependent function given by 

A IJK {K h, h) = £ £-g'« (-3^1 -^( kl + 2k 3 ) + h\ - kl kl\ 

perms * \ t / 

= £ irAS' K (-sf - + 2* 3 ) + i*{ - Ml) (69) 

perms V \ t / 

This is of order a slow-roll parameter, so (as we expect) the non-gaussianity will not be 
large when \e IJ \ <C 1. Eqs. (|68 )l - (|f)9jl represent the main result of this paper. 

The Feynman integrals that lead to (Q ! Q J Q K ) are insensitive to the field modes' 
behaviour deep inside the horizon and long after the modes have passed outside the 
horizon. The dominant contribution arises from fluctuations around the epoch of horizon 
crossing [£21] • For this reason, (|68|) represents the primordial three-point correlation 
among the scalars Q 1 a short time after horizon crossing. In particular, in making the 
estimate ([68 )1 - ([69 )1 . we have assumed that the three k-modes have roughly comparable 
wavenumbers, so that they cross the horizon at similar epochs. 

Eqs. (pjBjl - fJBT?)) do not by themselves represent the primordial non-gaussianity; 
instead, they only describe cubic interactions among the Q 1 at horizon crossing. In 
order to use Eqs. (ffi8' |) - ([B9"|) to calculate the primordial non-gaussianity in the curvature 
perturbation, it is necessary to use the general formalism outlined in Section El to 
assemble these initial non-gaussianities into the actual non-gaussianities that are in 
principle observable in the CMB. 

In the following sections, we consider two specific scenarios in which we can carry 
out this detailed assembly explicitly. We first verify that our results reduce to the 
well-known result of single-field inflation. We then proceed in Section El investigate the 
multiple-field analogue of the single-field power-law model known as assisted inflation 



4. Reduction to single-field case 

In this section we verify that Eq. (jo^j) reduces to the result of Maldacena [11 J for the 
case of inflation driven by a single scalar field. We work with the large-scale expression 
()23|). which requires an explicit expression for the number of e-folds, N, as a function of 
the field. For a single field <p with perturbation Q the integrated e-folding rate is given 
by (HZj): 

N = J Hdt = -^hi# (70) 

whenever e is approximately constant over the range of e-folds under consideration. 
Therefore, the derivatives of N, which are needed to determine (, are 

N =--^ N =-i5w + I^V (71) 
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Recalling that for single field models, the slow-roll parameters can be expressed as 



2 

H ) ' '' Hip ~ H 



* = 2(^V. ^ = -#7 = 2%, (72) 



we may write ( as 



/2l V 4 

As expected, to lowest-order this implies that the spectra are related by the well- 
known result [Hll E] 

(C(ki)c(ka)) = ^(g(k 1 )g(k 2 )), (74) 

with £ evaluated at the moment of horizon crossing. A similar expression holds for the 
three-point function. It follows from (|52*jl that the (CCC) an d (QQQ) correlators are 
related via 

(C(k 1 )C(k 2 )C(k 3 )) = -^^(Q(k 1 )Q(k 2 )Q(k 3 )) 

+ ^ (l " is") (^( k i)^( k 2)[Q*Q]( k 3)) + Perms, (75) 

where (as before) * denotes a convolution product. Eq. (|73j) implies that the connected 
part of the ( three-point function is related to the connected Q three-point function by 

(C(k 1 )c(k 2 )c(k 3 )) c = -^(g(k 1 )g(k 2 )g(k 3 )) c (76) 

where a subscript c denotes the connected part of a correlation function. However, 
the prescription of (f75]) means that we must mix some disconnected pieces with this 
connected component, which are described by the four-term correlators (QQQQ). These 
are no more than the superhorizon parts of the gauge transformation between the 
comoving and uniform curvature gauges, which were derived by Maldacena jTTj and 
found to have the form 

# 4 /1# 4 (p 1H 2 \ 1 



^(E^T{^m + iw)m m (77) 

Thus, the complete three-point function which we have calculated via this method is 
(C(k 1 )C(k 2 )C(k 3 )) = (2 7 r) 3 <K$>)^ TT i p -4, (78) 

, v 9 * iij zk i 

where A is given by 

* i>j y r * 7 i * i^j 

in exact agreement with Maldacena' s result. 
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5. Assisted inflation 

In this section we calculate the three-point correlation function for a multiple-field model 
of inflation known as assisted inflation |46J. In this scenario, a set of M scalar fields, each 
with an exponential potential, can conspire to drive a phase of accelerated expansion, 
even if the individual potential for each field is too steep to support inflation. The 
potentials are given by 



V I ( ¥ > I )=V exp (-y^V) 



(80) 




where Vo is a constant that determines the overall scale of the model, and p 1 represent 
coupling constants. + In the single field model, p 1 > 1 is required for inflation. The 
combined effective potential is given by V = J2i Vi- 
lli the spatially flat FRW cosmology, the system of fields ()8())1 admit a late-time 
attractor solution jHS] characterized by 

for some constants a 1 and any choice of reference field, which we have arbitrarily 
identified as if 1 . At the classical level, the existence of this attractor solution implies that 
the model is dynamically equivalent to that of a single scalar field with an exponential 
potential characterized by a coupling constant p = J2iP T ■ 

To proceed, we need to determine the number of e-folds in terms of the scalar field 
values. Following a standard calculation 0, we can express dN as 

dN = (82) 
<Pi 

for any field <£>j. When slow- roll is valid, we can also make use of the approximate field 
equations 3Hipj + dVi/dipi = 0, such that 

dN * dVjaV^ 1 " dlVc^ 1 = E dVjd^ d **- (83) 

The first and second derivatives of iV now follow immediately from Eq. (jHB^) : 
diV Vi 
dip j dVi/dipi 

Jrt-'-V'WfW,) ■ (84) 

Moreover, for exponential potentials it follows that 

Vj ./7 d"vj/dipj_ [T , sg) 



dVr/d^i V 2 ' dVi/dcp! \ p r 

+ The summation convention for indices J is suspended in this section. 
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Hence, the second derivative of N vanishes and the only remaining contribution to 
the total curvature perturbation (j2*Hj) is provided by the linear derivative term, which 
simplifies to 



< = -EVt '- (86) 

At the level of the two-point correlation function, this implies that 



2 



(C(ki)C(k 2 )> = J2 ^yM^V) = M^ki + k 2 )^ (87) 



(with an asterisk V indicating evaluation at horizon exit, as usual) and this is precisely 
the spectrum that would arise for a single-field model with an exponential potential 
characterized by p, in agreement with the classical observation that assisted inflation is 
dynamically equivalent to such a model [65J. 

On the other hand, we find that the three-point function is given by 



(c(k 1 )c(k 2 )c(k 2 )> = J2 v^^w^E^ttt^^* 1 '^'^- (88) 

I,J,L i Hi i 

In order to correctly estimate this amplitude, we need an expression for the slow-roll 
parameters. We take 

e 1 = (89) 
V 

Thus, tre /J = 1/p, which coincides with the slow-roll parameter in the classically 
equivalent single-field model. After some algebra it can be shown that this is equivalent 
to 

(c(koc(k 2 )c(k 2 )) = (2^^ ko=Ss ~p ( i £ k\k) -~5> 3 +^£ W) > (9°) 

i 1 % V * i>j i tyj J 

which is identical to the equivalent single-scalar model evaluated with Maldacena's 
single-field formula. 

This reduction to the single-field result illustrates a general feature of multiple-field 
models. In what amounts to the familiar Gram-Schmidt procedure, one can always 
make a rotation in field space so that one distinguished field (say Q l ) lies along the 
adiabatic direction, and the remaining M — 1 fields describe isocurvature fluctuations 
around it jSHl EH] ■ These isocurvature fluctuations do not contribute to the spectrum or 
bispectrum. Therefore, quite generally, one expects the resultant three-point function to 
be given by Maldacena's single-field result, constructed using the field-space trajectory 
of the adiabatic field. 

One should not understand this to mean that the formalism outlined in this paper 
is unnecessary. On the contrary, although the rotation in field space which produces an 
adiabatic field and M — 1 isocurvature fields can be performed explicitly in simple cases 
for a completely general multiple-field model, the most convenient way to account 
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for the effect of this rotation (at least on large scales) is to use exactly the non-linear 
extension [S^ of the Sasaki-Stewart 5N formalism, which was described in Section 12.21 
In this formalism it will always be necessary to 'prime' the superhorizon evolution using 
what amounts to an initial condition, describing the non-gaussianity produced at horizon 
crossing. Eq. (ffiHj) constitutes exactly the necessary initial condition. 

6. Conclusions 

In this paper, we have calculated the three-point function evaluated just after horizon 
exit that arises from a set of M scalar field fluctuations Q 1 coupled to gravity during 
inflation. The result is given by Eq. (|H8jl . where the momentum-dependent contribution 

JIJK 

is given in Eqs. ()68j) - (}oT?|) . This three-point function measures the intrinsic 
non-gaussianity produced at horizon exit in the coupled scalar-gravity system, and 
is important since it provides the initial condition needed for a calculation of the 
superhorizon evolution of the non-gaussianity in the inflationary density perturbation 
[TK] . Such superhorizon physics describes how the curvature perturbation (, which is 
used to calculate observable microwave background anisotropies via the Sachs-Wolfe 
effect, evolves outside the Hubble radius. Even though the initial non-gaussianity 
we have calculated, which is set at horizon exit, is typically small, the total non- 
gaussianity of ( after superhorizon evolution can be large. For example, this can occur 
in curvaton-like scenarios [HE] or where there is cross-correlation between the adiabatic 
and isocurvature modes (i(i . 

As a result, in some circumstances our expression could be used to verify that 
the non-gaussian effects generated at horizon crossing by microphysical processes are 
subdominant to the superhorizon, isocurvature-driven evolution of (, as described in 
|15j . However, more generally, both of these effects will be important and our result 
contributes in an essential way to the primordial non-gaussianity produced by inflation. 

It is important to note that (jfi8j) - (pTH|) describe the non-linearities among the Q 1 
at horizon crossing. The details of this were recently explored by Weinberg [S3], who 
showed that in a large class of multi-field models such integrals are dominated by the 
epoch of horizon-crossing, even beyond tree-level in perturbation theory. Provided one 
evaluates the slow-roll parameters when the k-modes in question crossed the horizon, 
Eq. (jnHJ) gives the relevant non-linearities. Of course, both the Q 1 and ( continue to 
evolve on superhorizon scales. For this reason one must assemble the (QQQ) correlators 
into (CCC) correlators to follow the time dependence correctly when outside the Hubble 
radius. The same conclusion was reached by Rigopoulos, Shellard & van Tent joTj . 

We do not address the question of what happens when one k-mode is squeezed, so 
that k\ <C k 2 , k s (say). In this case the mode corresponding to ki is pushed outside the 
horizon much earlier than the other two, since it is comparatively larger. In the single- 
field case, where perturbations "freeze in" outside the horizon owing to the constancy of 
(, the three-point function factorizes into the two-point function of k 2 and k 3 evaluated 
in a background which takes into account the back-reaction of the frozen ki perturbation 
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[TT] . In a multiple-field model, the k x mode will come to constitute part of the zero- 
momentum background when it is far outside the horizon, but since ( is no longer 
constant it is not so straightforward to provide a quantitative description of its effect. 
We leave this interesting question for future work. 

The level of non-gaussianity arising from microphysics is typically proportional to 
a slow-roll parameter, e 1 ^ 2 A IJK , and is therefore expected to be small when fluctuations 
are close to scale-invariance. Consequently, when isocurvature modes drive a very strong 
superhorizon evolution of (, and thereby source considerable non-gaussianity on very 
large scales, it seems likely that the generic situation will match that described by Lyth 
& Rodriguez ^H] , although the only model we are aware of where this can be verified in 
detail is the curvaton scenario. On the other hand, if the superhorizon contribution never 
becomes large, the non-gaussian effects we have considered will be important, in which 
case they can not be neglected. This is the situation that arises, for example, in the 
assisted inflationary scenario considered in Section^ which is the canonical model where 
the superhorizon piece is entirely absent. It is worth emphasizing that the three-point 
function can be calculated once the evolution of the unperturbed background cosmology, 
as parametrized by N^ 1 ), has been determined. Our result therefore complements that 
of and implies that whenever one has enough information about the background 
dynamics to calculate the momentum-independent piece of /nl, there will be always be 
enough information to calculate the momentum-dependent contribution as well. In this 
sense, no further work is required in order to calculate the full, momentum- dependent 
/nl- 

In general, the structure of the three-point function (|UH j) - (jp^j) is sensitive to the 
target space metric Q IJ . A similar dependence arises in the two-point function, which 
is proportional to Q IJ . Indeed, neglecting numerical coefficients, the form of (JoTJj) is 
in some sense inevitable. Specifically, it must arise from a cubic combination of the 
Q A, s which is a singlet with respect to the target space metric. Since the Q A can only 
appear contracted with the metric Qu, it follows that to leading-order in slow-roll, the 
only possible combination must have the form Gab<P a Q b GijQ i Q j ■ Such a combination 
leads to the index structure appearing in (joTJjl . there being no other target space vector 
at lowest-order in slow-roll which can be contracted with GabQ B to yield a singlet. 

It is natural to speculate that an analogous result will hold for the general n-th- 
order correlation function. When n is even (n = 2m) we expect the leading order 
slow-roll contribution to arise from a scalar quantity containing m contractions with 
Gij: 

/ . \ 2(m-l) 

( Q A Q B ---Q E Q F ) ~ f|J g^...g^l (91) 

2m copies " m copies 

where we have neglected momentum- dependent factors which may accompany each 
factor of the target space metric Q AB . When n is odd, so that n = 2m + 1, we anticipate 
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instead that 

/ • s 2(m-l) 

( Q A Q B Q C ---Q E Q F ) ~ I |J H- l ^ A g BC -..g EF \ . (92) 

2m + 1 copies m copies 

We have restricted our attention to the case where the target-space metric is 
flat and independent of the scalar field values, Q IJ = 5 IJ . Although this Ansatz 
covers a very wide class of multi-field inflationary models, it would be interesting to 
go beyond this approximation (see, for example, (Ml HE1 EH1 El])- Such an extension 
has been considered recently using a stochastic approach (THUS], although in principle, 
more general scenarios of this type could be investigated by employing the formalism 
developed in the present paper. For example, a specific point of interest concerns the 
curvature of the scalar field manifold. In general, a non-trivial metric Qu((f) will result 
in a curvature tensor of the form 

c\A A A F A F A /no \ 

" BCD = W BD,C — W BC,D + W BD^ FC ~ W BC^ FD, (."<jj 

where uj A bc is the Levi-Civita connexion compatible with Qjj. This tensor is identically 
zero for the models considered in this paper, but is expected to arise in the expressions 
for the spectral tilt and other cosmological observables 03 m more general classes 
of models. It is possible that the presence of curvature terms of this form may be 
sufficiently important to invalidate the flat target-space analysis we have adopted. In 
this larger primordial non-gaussianity could be generated. However, since the 

curvature tensor VL A bcd can not be too large if slow-roll is to be respected |^H] . we 
anticipate that this will not be the case in general. 
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